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FOREWORD 


The  work  described  in  this  report  was  done  during  the  period  from 
September  1972  through  June  1973  under  Contract  Number  F08635-71-C-0090 
with  the  Air  Force  Armament  Laboratory,  Eglin  Air  Force  Base,  Florida. 

The  program  managers  were  Major  William  A.  Miller  and  Captain  Visi 
Arajs  (DLJA).  This  work  was  part  of  a  follow-on  effort  of  the  work  re¬ 
ported  as  AFATL-TR-71-109.  Cross-Flow  Corrected  Axisymemtric  Solution  for 
Multiple  Body  Interferences,  dated  August  1971. 

The  Computer  program  for  this  work  was  written  by  Charles  J.  Smith 
and  Grady  H.  Saunders,  Jr.,  Graduate  Assistants,  and  Dr.  M.  A.  Cutchins, 
Associate  Professor,  Aerospace  Engineering,  Auburn  University.  All 
experimental  work  under  this  project  was  accomplished  by  the  Department 
of  Aerospace  Engineering,  Auburn  University,  Auburn,  Alabama. 
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ABSTRACT 


The  mutual  aerodynamic  interference  problem  for  two  axisymmetric 
bodies  has  been  analyzed  using  the  image  system  technique.  In  order  to 
facilitate  this  analysis ,  it  has  been  assumed  that  small  perturbation 
solutions  Eire  vEilid.  It  is  further  assumed  that  the  exterae!  stores  are 
slender  bodies  and  that  the  mutusil  interference  cem  be  analyzed  by  fi  rBt 
assuming  a  cross-flow  solution.  T~a  image  system  in  the  cross-flow 
plane  consists  of  source-sink  pairs  appropriately  located  by  using  the 
Milne-Thomson  circle  theorem.  The  actuELL  three-dimensionELl  source-sink 
pairs  are  displaced  from  the  body  axis  according  to  the  cross-flow  image 
system.  Their  strengths  are  then  determined  by  the  Rankine  method. 

Good  agreement  has  been  found  between  the  theoretics!  and  experimental 
results . 
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SECTION  I 


INTRODUCTION 


This  work  is  part  of  a  continuing  effort  to  analytically  examine  the 
flow  field  beneath  an  aircraft  with  external  stores  in  an  effort  to 
produce  equations  for  computing  store  separation  trajectories  which  include 
the  effects  of  mutual  aerodynamic  interference.  The  application  o£. the 
cross-flow  solution  is  an  alternate  approach  to  the  time-consuming  and  very 
complex  method  of  solving  the  interference  problem  by  using  a  prohibitively 
large  number  of  singular  solutions,  such  as  that  obtained  by  Smith  and 
Pierce  (l)  for  a  surface  source  distribution  on  all  the  bodies.  Thus, 
the  most  obvious  advantage  of  the  cross-flow  correction  method  is  its 
relatively  simple  application  to  the  case  of  interfering  multiple  dissimilar 
bodies . 

The  previous  analytical  solution  proposed  by  Martin  (2)  in  which  the 
mutual  aerodynamic  interference  problem  was  considered  by  using  a  cross- 
flow  corrected  axisymmetric  solution,  was  re-evaluated.  In  that  work  an 
approximate  line  source  distribution  was  assumed  as  the  solution  for  the 
isolated  axisymmetric  body.  This  solution  was  then  used  to  generate  the 
appropriate  image  system  as  suggested  by  the  two-dimensional  cross-flow 
solution  outlined  in  Reference  3.  In  order  to  obtain  a  better  correlation 
with  experiment  and  to  generate  a  consistent  technique  that  could  be  used 
for  any  axisymmetric  body  shape,  this  more  general  approach  has  been 
developed. 

The  calculation  of  the  flow  field  for  small  disturbances  requires 
a  solution  to  the  Prandtl-Glauert  equation  subject  to  the  boundary  condi¬ 
tions  of  no  fluid  flow  into  the  body  and  vanishing  perturbations  at  large 
distances  from  the  body.  Since  this  equation  is  linear,  superposition  of 
elementary  solutions  will  allow  the  build-up  of  a  complex  flow  which  satis¬ 
fies  the  boundary  conditions.  The  widely  used  vortex-lattice  and  line 
source  solutions  for  planar  and  axisymmetric  flows  represent  two  common 
examples  of  the  superposition  of  simple  solutions  to  form  more  complex 
flow  fields .  However ,  when  systems ,  such  as  the  flow  fields  for  planar 
wings  and  axisymmetric  bodies,  are  combined,  it  is  not  possible  to  simply 
superimpose  the  results  for  the  separate  systems .  The  interference  of  the 
flow  fields  requires  that  adjustments  be  made  to  the  elementary  source 
and/or  vortex  distributions  in  order  to  satisfy  the  body  boundary  conditions . 
The  cross-flow  correction  discussed  herein  is  an  example  of  this  type  of 
adjustment  for  two  interfering  axisymmetric  bodies . 

This  solution  is  based  on  a  two-dimensional  image  system  which  utilizes 
the  Rankine  method  (U)  of  superposition  of  an  axial  source-sink  distribu¬ 
tion  to  form  a  closed  body.  In  the  cross-flow  solution,  however,  the 
sources  and  sinks  are  displaced  from  the  axes  of  the  bodies  in  such  a 
manner  as  to  account  for  the  interference  effects  and  to  satisfy  the  bound¬ 
ary  conditions.  The  derivation  of  the  cross-flow  image  system  (3) 
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utilizes  the  Milne-Thomson  circle  theorem  (5)  and  is  based  on  the  premise 
that  the  boundary  condition  for  two-dimensional  interference  flow  past 
two  circular  cylinders  can  be  satisfied  by  displacing  the  sources  in  each 
body  towards  each  other  by  a  specified  distance. 

This  solution  differs  from  that  of  Reference  2  in  that  instead  of 
assuming  that  the  strengths  of  the  sources  in  the  image  system  are  known, 
as  in  the  previous  work,  these  values  are  determined  so  that  when  added 
to  the  predetermined  axial  source-sink  distribution,  the  body  boundary 
conditions  are  satisfied. 

The  new  technique  was  applied  to  calculate  the  interference  flow 
field  for  two  M-117  bombs.  Results  of  the  theoretical  calculations  were 
compared  with  the  experimental  values  as  obtained  by  tests  performed  in 
the  low-speed  wind  tunnel  at  the  contractor  facility. 

The  computer  program  used  to  compute  the  theoretical  values  is  given 
as  Appendix  I. 
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SECTION  II 


THEORETICAL  CONSIDERATIONS 


The  method  of  solution  for  analysis  of  the  two-body  aerodynamic  inter¬ 
ference  problem,  as  developed  by  Martin  in  Reference  3,  has  been  investiga¬ 
ted,  and  a  new  technique  for  improving  this  method  has  been  developed.  The 
basis  of  Martin's  method  was  the  combination  of  the  single  body  axisymmetric 
solution  and  the  cross-flow  image  system  [3]  to  approximate  the  interference 
flow  field  for  two  similar  bodies  of  revolution  in  a  uniform  flow.  It  was 
found,  however,  that  this  technique  does  not  satisfy  the  boundary  condi¬ 
tions  for  the  three-dimensional  problem  [6],  In  order  to  correct  the  body 
boundary  conditions ,  the  axisymmetric  solution  was  combined  with  a  cross- 
flow  solution  consisting  of  a  source  distribution  displaced  from  the  body 
center  line  and  a  sink  distribution  of  equal  magnitude  along  the  body 
axis  as  suggested  by  the  two-dimensional  cross-flow  solution. 

The  usual  assumptions  of  small  disturbances  and  inviscid  flow  are 
assumed.  Further,  it  is  assumed  that  the  flow  is  incompressible  so  that 
the  familiar  Laplace  equation  is  the  governing  equation  and  the  well-known 
elementary  solutions,  source-sink ,  doublet,  etc.,  flows  areavailable. 

Since  the  body  boundary  condition  for  the  two-body  problem  cannot  be 
satisfied  by  simply  superimposing  the  solutions  for  isolated  bodies , 

Martin  [2]  developed  the  cross-flow  image  system  as  a  possible  correction. 

In  his  original  work,  Martin  proposed  a  solution  in  terms  of  a  continuous 
distributed  line  source  which  was  displaced  according  to  the  two-dimensional 
image  system.  Generally,  it  is  not  possible  to  find  a  continuous  line 
source  for  generating  the  body,  and  the  solution  is  approximated  by  a  number 
of  finite  line  sources.  This  approach  leads  to  geometric  complexities 
when  the  image  system  is  applied,  because  the  line  sources  must  be  at 
angles  corresponding  to  the  slope  of  the  displacement  curve  (Figure  l). 

For  this  reason,  the  line  sources  were  replaced  with  a  number  of  three- 
dimensional  point  sources  and  sinks.  This  also  has  the  advantage  of  re¬ 
quiring  half  as  many  calculations  as  the  solution  for  the  same  number  of 
line  sources,  because  the  integrated  effect  of  a  line  source  at  a  field 
point  is  a  function  of  the  distance  from  each  end  of  the  source  to  the 
point,  whereas  the  effect  of  a  point  source  depends  only  on  the  distance 
between  the  points . 

Any  particular  solution  of  Laplace's  equation. 


iii  + 

ix2 


ill  +  lii 

3y2  3z2 


-  0 


(1) 


-3- 


-4- 


Figure  1.  Image  System  Representation  o:'  the  Tvo-Body  Case 


which  is  representative  of  an  actual  flow  must  satisfy  the  boundary  con¬ 
ditions  at  infinity  and  at  the  surfaces  of  any  bodies  which  are  present 
in  the  field.  I  physically  meaningful  boundary  condition  is  that  the 
perturbation  velocity  components  tend  towards  zero  at  large  distances  from 
the  disturbances.  This  is  known  to  be  the  case  for  source  flow,  as  will 
be  evident  from  the  velocity  derived  from  the  potential  function.  The 
boundary  condition  along  the  surface  of  an  impermeable  body  must  express 
the  condition  of  no  fluid  flow  into  the  body.  This  condition  is  satis¬ 
fied  when  the  component  of  velocity  normal  to  the  body  is  zero.  Equivalent¬ 
ly,  the  boundary  condition  is  satisfied  when  the  velocity  at  the  fluid- 
solid  interface  is  tangent  to  the  body.  For  an  axisymmetric  body,  the 
boundary  condition  may  be  expressed  mathematically  as 


dr 

dx 


(2) 


where  dr/dx  is  the  slope  of  the  body  in  any  meridional  plane,  Uj.  and 

are  the  perturbation  velocity  components  in  the  radial  and  social  directions, 

and  U„  is  the  magnitude  of  the  free-stream  velocity. 

For  a  single  body  at  zero  angle-of-attack  the  angular  velocity  com¬ 
ponent,  ug,  is  identically  zero.  In  the  two-body  case,  u©  is  not  generally 
zero,  but  it  does  satisfy  the  tangential  flow  requirement  at  all  points 
on  an  axisymmetric  body.  Therefore,  equation  (2)  is  the  correct  expression 
for  the  body  boundary  condition  for  the  two-body  case  as  well  as  for  the 
isolated  body. 


-5- 


SECTION  III 


SINGLE  BODY  SOLUTION 


Assuming  that  any  axi symmetric  body  can  he  represented  by  sources 
distributed  along  its  center  line,  the  total  velocity  potential  becomes 
the  sum  of  these  sources  plus  a  constant  and  can  be  expressed  as 


N  m 

<f>,  =  -  2  t—  - —  +  4i  *  (3) 

J  i  1471  Pij 

where  is  the  distance  from  the  ith  source  to  any  point  J  in  the  flow 
field  ana  <f>’  is  a  constant  and  represents  the  free-stream  velocity  potential. 


Differentiating  the  summation  term  in  equation  (3)  yields  the  per¬ 
turbation  velocity  components  which,  in  cylindrical  coordinates  r,  0,  x, 
can  be  expressed  as 


N 

u  *  -  I 
XJ  i 

mi 

77 

—  (—  ) 

9XJ  PiJ 

(h) 

N 

u  =  -  E 

rJ  i 

mi 

77 

If  (  r- 1 

8rj  pu 

(5) 

N 

X  '  - 

mi 
hTtr . 

L.  (_L_) 

90,  V/ 

(6) 
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Due  to  symmetry,  u^  is  zero  for  the  case  of  a  single  body  at  zero 
angle-of-attack.  As  can  be  seen  from  Figure  2,  p..  is  not  a  function  of 
9  and  is  given  by  J 


P 


ij 


(7) 


If  a  new  variable  z^ , 
source  strength,  i.e., 


defined  such  that  it 


is  proportional  to  the 
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Figure  2.  Geometric  Relations  for  the  Single  Body  Case. 


(8) 


n. 


z,  = 


*i 

Wu 


is  introduced  and  equations  CM  and  (5)  are  used,  the  velocity  components 
can  be  written  in  the  forms , 


(9) 


N 

+  E  z 
i 


ij 


(10) 


Equations  (9)  and  (10)  can  be  applied  at  any  non-singular  point  in  the 
flow  field  to  determine  the  velocity  ratios.  Referring  to  equation  (2), 
it  can  be  seen  that  the  body  boundary  condition  in  terms  of  the  velocity 
ratios  is  specified  by 
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ii 
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(11) 


By  multiplying  both  sides  of  equation  (ll)  by  the  denominator  on  the 
right-hand  side  of  that  equation  and  simplifying  the  resulting  equation, 
one  may  obtain  the  result, 


N 

E 

i 


z. 

i 


<yV 


dr 


dx. 


+  r. 


p! 


ij 


( J=1 ,2 ,3, . . . ,N) 


(12) 


Equation  (12)  forms  a  set  of  n  simultaneous  linear  equations  for 
determining  the  required  source  strengths  in  terms  of  the  known  geometry. 
If  the  coefficients  of  z.  in  equation  (12)  are  defined  to  be  c.  .,  the 
equation  may  be  written  In  the  form  J 


N  dr, 

E  zicij  =  d x  (J=1,2,3,...,N)  (13) 

i  J 
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Equation  ( 13)  may  be  written  in  the  equivalent  matrix  form. 


[C](Z)  =  (dr/dx) 


(lM 


Solving  for  (z)  by  matrix  inversion,  one  obtains  the  result 

(Z)  =  [cHUr/dx)  (15) 

The  array  of  strengths  as  calculated  by  this  method  satisfies  the  body 
boundary  condition  at  the  control  points  for  any  isolated  axisymmetric 
body  at  zero  angle-of-attack. 

The  predicted  pressure  distribution  for  an  M-117  bomb  is  compared 
with  experiment  in  Figure  3  and  shows  good  agreement. 


-9- 


SECTION  IV 


SOLUTION  FOR  TWO  SIMILAR  BODIES 


Consider  the  case  of  two  identical  axLsymmetric  bodies  aligned  such 
that  one  is  directly  above  or  to  the  side  of  the  other.  According  to 
Reference  2,  the  image  system  for  this  case  is  formed  by  displacing  each 
source  towards  its  corresponding  source  in  the  other  body,  as  in  Figure 
U.  The  distances  which  the  sources  are  displaced  are  derived  from 
repeated  applications  of  the  Milne-Thomson  circle  theorem  [3]  and  are 
given  by 

6±  =  i  {l-[l-C2ri/yo)2],S}  (l6) 

where  yQis  the  distance  between  body  center  lines  and  r^  is  the  radius  of 
the  body  at  the  ith  source.  Since  the  bodies  are  identical,  the  sources 
at  corresponding  axial  locations  in  each  body  are  of  the  same  strength. 

In  the  two-dimensional  case,  the  image  system  consists  of  the  dis¬ 
placed  sources  with  corresponding  sinks  of  equal  strength  at  the  center 
line.  Since  the  strengths  of  the  image  pairs  at  any  cross-section  are 
the  same  as  that  of  the  original  source,  from  the  single  body  solution, 
the  sink  cancels  this  source.  Hence,  the  two-dimens ional  system  yields 
only  the  displaced  line  of  sources  as  discussed  above.  In  the  three- 
dimensional  case,  however,  this  simple  system  does  not  produce  a  satis¬ 
factory  result. 

In  the  three-dimensional  case,  there  is  a  significant  component  of 
flow  in  the  axial  direction  in  addition  to  the  radial  flow.  It  is 
reasoned,  however,  that  an  equivalent  line  source  can  be  assumed  at  the 
body  center  line  which  contributes  only  to  the  radial  flow.  This  is 
purely  a  fictitious  flow  but  suffices  to  justify  an  equivalent  two- 
dimensional  flow  field  which  would  generate  the  sou-ce-sink  pair  described 
in  the  image  system  above.  Note  that  the  source  displacement  distance, 

6,  depends  only  on  the  geometric  relations  and  is  independent  of  the  two- 
dimensional  source  strength. 

It  follows  that  an  appropriate  system  for  the  three-dimensional  case 
can  be  generated  by  retaining  the  sources  along  the  center  line  of  the 
two  bodies  (isolated  solution)  and  adding  an  image  system  consisting  of 
source-sink  pairs  in  each  body  which  are  located  according  to  the  two- 
dimensional  displacement  distances. 
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The  total  velocity  potential  for  this  two-body  system  can  now  be 
expressed  as  before  and,  for  a  body  formed  by  N  sources,  the  potential  at 
a  point  is 


N 

-  1 
i=l 


ItirV' 


ij 


(17) 


where  m^  refers  to  the  single  body  source  strength  and  hl[  is  the  strength 
of  the  added  source  and  sink.  p.  ,  pp  ,  p'  ,  and  p'  are  as 

ij  ij  ij  ij 

illustrated  in  Figure  4  and  may  be  expressed  as 


P1  =  [(<5i"rj  sin  9)2  +  (xj"Xi)2  +  (rJ  cos  e)2]*5  (18) 

i  J 

p2  =  ^yo_(Si_rJ  sin  e^2  +  (xj”xi )2  +  (rj  cos  ^)2]Js  (19) 

^  J 

p'  =  [(x  -x  )2  +  r  2]**  (20) 

xij  J  1  3 

Pp  =  [(y  -r.  sin  e)2  +  (x,-x, )2  +  (r  cos  e)2]5*  (2l) 

2, ,  °  J  J  i  J 


Using  equations  (4),  (5)  and  (6)  for  the  velocity  components  and 
equation  (8)  for  zis  one  may  obtain  the  velocity  ratios 
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)] 


Ur  N 


u 


8.  ,  N 

=  -  —  z  z  d!,  +  z!  (d,,  -  d'  )] 
U»  rJ  i=l  1  lJ  1  iJ  iJ 


_  3  t  1  .  1 _ v 

"  3x  p  p 

3  XiJ  ij 


(22) 


(23) 

(24) 


(25) 
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—  (- 
30 


ij 


(26) 

(27) 


The  terms  a'^,  b'^.  811(1  8X6  forne<i  from  equations  (25) »  (26),  and 

(27)  with  the  values  for  replaced  with  the  primed  quantities  (see 

Figure  1+  and  equations  (l8)  to  (2l)).  The  unprimed  values  of  z  correspond 
to  the  single  body  source  strengths ,  and  the  primed  values  represent  the 
image  system. 


The  source  distribution  to  correct  for  the  interference  is  developed 
by  application  of  the  requirement  of  tangential  flow  at  the  body.  Using 
equations  (22)  and  (23)  for  the  velocity  components  and  equation  (2)  for 
the  boundary  condition,  one  may  obtain  the  result 


dr 

dx 


=  l  (z,( 


dr 


J 


i  ij  dx. 


-b|j)  *  *it(. 


dr 


ij  dx, 


-v-< 


dr 


ij  dx 


J 


-b'j)]}  (28) 


Equation  (28)  represents  a  set  of  N  simultaneous  equations  for  determining 
the  source  distribution  required  to  account  for  the  interference  between  two 
similar  bodies.  By  defining  the  coefficients  of  z.^  and  z1.  as  e^  and  f^  , 
respectively,  equation  (28)  can  be  written  in  the  equivalent  matrix  form 

(dr/dx)  =  [e](z)  +  [f](z')  (29) 

Solving  this  equation  for  ( z ' )  one  obtains 

(z1)  =  [f]”1  {(dr/dx)  -  [e](z)>  (30) 

The  source  distribution  given  by  equation  (30)  will  satisfy  the  body 
boundary  conditions  at  N  control  points  and  give  good  agreement  at  all 
other  points  around  the  body.  It  was  found  that  the  best  location  for 
the  control  points  lies  along  the  intersection  of  one  body  with  the  plane 
which  passes  through  the  center  lines  of  both  bodies;  i.e.  ,  from  Figure  4, 
the  90-degree  meridional  plane.  The  predicted  slopes  at  other  meridional 
planes  are  shown  on  Figure  5  for  the  M-117  bomb  body  shape . 

The  predicted  pressure  distribution  and  experimental  points  are  shown 
in  Figures  6  to  11  for  meridional  angles  between  90°  and  -90°.  These 
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distributions  were  calculated  using  a  30-point  solution  for  the  axisym- 
metric  body  and  for  the  cross-flow  correction. 


90°  (Actual  Body  Slope) 


Figure  5.  Comparison  of  Body  Slope  with  Computed  Streamlines 
at  Meridional  Angles,  0  =  -  90°,  0°,  90°. 
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Lgure  6.  Comparison  of  Theoretical  and  Experimental  Pressure  Distribution 
for  the  Two-Body  Case,  y  =  1.0U2  Dia. ,  0  =  90°. 
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Figure  7-  Comparison  of  Theoretical  and  Experimental  Pressure  Distribution 
for  the  Two-Body  Case,  yQ  =  1.0^2  Dia. ,  0  =  60°. 
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Figure  10.  Comparison  of  Theoretical  and  Experimental  Pressure  Distribution 
for  the  Tvo-Body  Case,  y  =  1.0^2  dia. ,  0  =  -  30°. 
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Figure  11.  Comparison  of  Theoretical  and  Experimental  Pressure  Distribution 
for  the  Two-Body  Case,  y  =  1.0^2  Dia. ,  0  =  -  90°. 


SECTION  V 


SOLUTION  FOR  TWO  DISSIMILAR  BODIES 


The  two-body  solution  of  the  previous  section  can  be  extended  to 
include  the  case  of  dissimilar  axially  misaligned  bodies.  Again,  and  in 
the  solution  for  two  similar  bodies,  it  will  be  assumed  that  the  isolated 
body  solution  has  an  equivalent  two-dimensional  source  distribution  which 
causes  the  radial  flow  outward  from  the  body  centerline.  Also,  it  should 
be  noted  that  this  radial  flow  occurs  both  forward  and  aft  of  the  body 
as  well  as  along  the  body.  Thus,  in  order  to  generate  the  cross-flow 
image  system,  it  is  assumed,  conceptually,  that  the  flow  field  forward 
and  aft  of  an  isolated  body  is  spread  by  equivalent  sources  distributed 
along  the  extension  of  the  body  center  line.  The  image  system  in  body 
2  at  cross- flow  station  a  (see  Figure  12)  can  be  generated  using  the 
Milne-Thomson  circle  theorem  for  source  flow  over  a  cylinder  and  is  illu¬ 
strated  as  Figrre  13(a)  where  is  given  by 


y0 


The  complete  image  system  at  station  b  can  be  generated  by  repeated 
application  of  the  Milne-Thomson  circle  theorem  and  consists  of  a  line  of 
sources  and  sinks  within  both  bodies.  The  complete  development  is  given 
in  Reference  3.  It  has  been  found,  however,  that  one  source-sink  for 
each  body  with  the  source  displacement  distances  given  by 


2 
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Figure  12.  Geometric  Relations  for  Two  Misaligned  Bodies. 
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and 


yields  good  results  as  will  be  shown  later.  Also  note  that  equations  (3?) 
and  (3^)  reduce  to  the  truncated  series  representing  equation  (l6)  for  6 
in  the  solution  for  two  similar  bodies  where  r^  =  r^. 

Because  of  the  asymmetry  of  this  case,  the  strength  and  the  dis¬ 
placements  of  the  source-sink  pairs  are  different  for  the  two  bodies. 

Thus,  if  there  are  N  image  pairs  in  each  body  (corresponds  to  N  sources 
for  the  isolated  body  solution),  then  there  will  be  a  set  of  2N  simultan¬ 
eous  equations  to  be  solved. 

The  velocity  potential  function  for  this  case  can  be  expressed  as 


N 

-  Z 

i=l 


N 

I 

i=l 


N 

Z 

i=l 


( 


N 

Z 

i=l 


ij 


ij 


(3b) 


where  the  subscripts  1  and  2  correspond  to  bodies  1  and  2,  respectively, 
as  shown  on  Figure  12.  The  other  terms  are  consistent  with  the  previous 
notation  and  as  defined  by  Figure  12. 

Performing  the  appropriate  differentiations  on  equation  (34)  in 
order  to  obtain  the  velocity  components  and  using  equation  (2)  for  the 
boundary  condition  generates  the  set  of  2N  simultaneous  equations  which 
can  be  written  in  matrix  form  as  equation  (30).  The  z'  values  in  this 
case  are  N  values  of  z'^  and  N  values  of  z'2  corresponding  to  the  image 
source-sink  pairs  in  bodies  1  and  2,  respectively. 

The  N  control  points  on  body  1  were  taken  along  the  top  of  the  body 
at  the  intersection  of  the  meridional  plane  passing  through  both  body 
center  lines.  The  other  N  control  points  were  taken  at  the  corresponding 
lower  surface  of  body  2. 
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The  predicted  pressure  distribution  and  experimental  results  for 
body  1  with  body  1  displaced  two  diameters  downstream  from  body  ?  are 
shown  on  Figure  lU. 
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Comparison  of  Theoretical  ai 
Distribution  for  the  Solui  i( 
=  °.0  Dia.  and  y  =  1  .04; 


x 


SECTION  VI 


CONCLUSIONS 


This  method  led  to  good  agremeent  between  the  experimental  and 
theoretical  pressure  distributions.  Also,  the  slope  of  the  body  stream¬ 
line  calculated  using  this  method  and  the  slope  of  the  actual  body  were 
in  good  agreement.  Thus,  this  method  can  be  considered  a  valid  solution 
for  the  mutual  aerodynamic  interference  problem  for  two  axisymmetric 
bodies  at  zero  angle-of-attack. 
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appendix  I 


FORTRAN  COMPUTER  PROGRAM  FOR 
2-BODY  INTERFERENCE 


The  program  which  follows  is  the  FORTRAN  computer  program  used  to 
calculate  the  flow  field  for  two  parallel  bodies  in  close  proximity  co 
each  other  by  the  source-sink  image  pairs  approach.  The  program  was 
written  for  the  IBM  360  system  and  compiled  on  the  FORTRAN  G  Compiler. 
All  calculations  were  performed  in  single  precision  arithmetic. 
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ooo  ooo  n  ooo 


fWO  SIMILAR  BODIES 


DIMENSION  PAD (31 ) ,X(3I ) , SLOPE ( 31) , XA( 30>,XB( 30  ) , DRDXA <  30 ) , 
-Z  A ( 3 J )  *  DPCXB ( 30),DELA(30)»D£LB(30) 

DIMENSION  R  A  (  3  0  )  »  RB(30),  ZB(30) 

COMMON  /TWC  /  XA, X0,RA, DRDXA, CRDXB , DE L A , D t LR , 7  A 
EQUIVALENCE  ( R A ( I ) v  RB (  l ) , R AD  (  1  )  ) ,  ( L A ( 1 ) , L B (  1 )  ) 


NS  =  30 
DX  =  0 .20 
X0=0.0 

CALL  GtCM(RAC,X, DX, NS, SLOPE  ) 

400  CONTINUE 

I F (  XO  .NE.  0.0  )  GO  TC  18 

SINGLE  BODY  SOURCE  STRENGTHS 

Y0  =  1 0000. 0 
GO  TO  17 

TWO-BODY  SOLUTION 

18  CONTINUE 
Y0  =  1 .042 

17  DC  2C0  1=1, NS 
X A (  I  )=X(I  )  +  X0 
XB ( I  )  =X ( I ) 

DRDXA ( I )  =  SLCP£ (  I  ) 

DR CX B ( I )=-SLOPE(  I  ) 

200  CGNTINUt 

IF (YO.Gt. 10000.0)  GO  TC  16 

DC  150  1=1, NS 

M=  X  A  (  I  )/DX 

N= ( XB ( I )-XC)/DX 

IF(M.GT.NS)  GO  TO  156 

0ELA{ I )=RA(  I )**2/(Y0-( RB(M)**2)/Y0) 

156  CONTINUE 

IF ( X A ( I ).Gt.XB(NS)  )  DELAll )=RA(I)**2/Y0 

IF(N.Lt.O)  GC  TO  155 

OELB ( I ) =RB ( I )**2/(Y0-(RA(N)**2)/Y0) 

155  CONTINUE 

I F ( XB ( I ).LT.XO)  DELB (  I  ) =RB ( I  >**2/Y0 
150  CONTINUE 

WR  I  TE ( 6, 1 1C )  XO*  YO 

110  FORMAT ( 1 H 1 , •  TWO-BODY  SOLUTION'/'  XO  =  ',F10.5,*  YO  =  ', 
-f 10. 5//  '  I  PELA(I)  DELB ( I ) • // ) 
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WRITL<6,111M  I , DEL  AC I) ,CELB(  I  ), 1  =  1, NS) 
111  FORMATC IX, 12.2F12.5/) 

16  CALL  CROSS 2 ( NS , XO , YO ) 

IFCYO.GE. 10000.0)  GO  TO  18 
I F (  XO  .NE.0.0  )  GO  TO  300 
X0=2.0 
GO  TO  ACO 
300  STOP 
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noon 


C 


SUBROUTINE  CRO SS2 ( N S, XO , YO ) 


DIMENSION  ZCO(60,60),  EI60),  DELAI 30  ) , DELB  (  30 > ,  X ( 30  >  * 

R  1 30 ) , X A I  30 ) ,  XB(  30 ) , DRDX A I  30 ) , DRDXB  I  30 ) ,ZA(30) 

-  ,UX(30),UR(30),BCI30),CP(30> 

DIMENSION  PC  1 30, 30 ) »  KAI30),  RB 1 30 ) , Z B I  30 ) » Z AP l 30 ) , Z BP  I  30  ) 
DIMENSION  0 T ( 30 ) 

COMMON  /TWO/  XA,XB,RA,DRDXA, DRDXB, DELA, DELB, ZA 
tQUlVALENCfc  (PCI  1,1) ,ZCO( 1,1)1, 

(KA(I)tRBI  I)),(ZA(I),ZB(  l)), 

(ZAPm,E(l)),  IZBPI1),  r.(  31)  ) 

S I ND  f  PPP ) =  SINIPPP/57. 29578) 

CGSQ(PPP)=  COSIPPP/57. 29578) 


EVALUATION  ALONG  BODY  A 

DO  301  J=1,NS 
X  (  J )  =XA ( J ) 

R  (  J ) =RA ( J  ) 

E  (  J ) =  DKDX  A ( J  ) 

301  CONTINUE 

DO  100  J*  1  ,NS 
00  100  1=1, NS 

P  A  1  =  SORT!  I  X  I J ) -X A ( I ) )**2+RIJ)**2) 

11  AX=(X{J)-XA(I)  )/PAl**3 

12  AR=R< J)/PA1**3 

IF (YO.LT. ICOOO.O)  GO  TO  10 
PCI J, I )=AR-AX*DRDXA( J) 

GO  TO  100 
C 

10  P A2=  SORT!  (X(J)-XA(I) )*+2+(DELA( I )-R I J ) )**2 ) 

PB 1 =  SORT! I  X ( J )-XB  (  I ) )**2+(R(J)-Y0)**2) 

PB2=  SORT ( (XI J )-XB(  I ) ) **2+ IR  I  J ) +DELB  I  I )-Y0)**2) 
C 

BX=( X< J)-XA(I) )/PA2**3-iX( J)-XA( I  )  )/PAl**3 
CX  =  ( X ( J )-XB ( I )  )/PBl**3 

DX=(X( J )-XB( I) )/PB2**3-(X( J)-XB(I ) )/PBl**3 

C 

BK=(R(J) -DELAI I ) )/PA2**3-R(J)/PAl**3 
CR=(  I J  >-Y0  )/PBl**3 

DR  =  (KI J  >*DELB(  I )-YO)/PB2**3- IRI J)-Y0)/PBl**3 
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C 


aa=ar-ax*drcxa ( J) 


o  o  o  o 


CA=CR-CX*DRCXA< J ) 

C 

ZCQI J , I  )=BR-BX*DRDXA( J ) 

ZCG( J, I+NS)=DR-DX*DRDXA( J) 

EIJ)  =  fc I J )-Z A (  I  ) *A A-ZB (  I  )*CA 

C 

ICO  CONTI NUfc 
C 

IFIYO.LT. 10000.0)  GO  TO  14 
CALL  SIMQ  (PC,DRDXA,NS,KS) 

DO  15  I =1 ,NS 
15  Z A { I )=DRDXA( I ) 

KtTURN 

C 

14  CONTINUE 


EVALUATION  ALONG  BODY  B 

DO  302  J=  1  » N S 
XI  J)=Xl5(J) 

R ( J ) =YO-RB ( J ) 

E( J+NS)-DRCXB( J) 

302  CONTINUE 
C 

DO  200  J-  1 1 NS 
DO  2C0  1=1, NS 
C 

PA  1=  SORT! (XIJI-XAI  !)) **2+R ( J ) **2 ) 

PA2=  SORT! (XI JJ-XAI  I  ) ) **2+ 1  DEL A ( I )  —  R ( J I i**2 ) 

PB 1=  SORT! (X (J)-XB(  I  ) )**2HR(  J)-Y0)**2) 

P62=  SURTI  (XUI-Xfil  I  ) )**2+(R( JI+DELBI I)-Y0)**2) 

C 

A  X  * ( X ( J  )-XA(  I  ))/PAi**3 

HX=(X(J)-XA(I  ) )/PA2**3-(X(J)-XAI I J  ) /PA  1**3 
CX  =  (X(J)-XB(  I)  ) /PB 1 **3 

l)X=  I  X  (  J  )  -XB  I  I  )  )/PB2**3-(X(  JI-XBI  I  )  >/PBl**3 

C 

AR=R ( J ) /PA  l  **3 

PR=(R(  J  )-CELA(  m/PA2**3-R(  J) /PA  1**3 
CR*(RI J ) — Y  0  )  /PB 1**3 

DR  =  IR« J I  +  DELBl 1 1 — VO  )/PB2**3-lR(J)-Y0)/PBl**3 

C 

C 

AB=AR-AX*DRCXBI J ) 

CB=CR-CX*DRDXB< J) 

C 

ZCO( J ♦ N S , I  )  =8R-BX*DRDXB( J) 
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O  O  O  O  o  o  o 


Z  CC  (  J*NS,  I+NS)=DR-DX*DRDXB( J) 
fc<  J  +  NS) =  E  l  J+NSJ-ZA ( I ) *AB-ZB ( I ) *CB 

C 

200  CONTINUE 
NS2=2*NS 
C 

CALL  SIMQ  (  ZCO  1 1 » NS2  »  K  S ) 

IHQS.tC.  1.0)  GO  TO  2200 
GC  TO  300 
2200  WR  I  TE ( 6  »  40C ) 

400  FORMAT! IH1, 'SOLUTION  FOR  SOURCE  STPENGHTS  IS  SINGULAR*,/) 
RETURN 

300  WR I T  E ( 6 , 500 ) 

300  FORMAT!  1H1,//,  IX, 'NON-DIMENSIONAL  SOURCE  STRENGHTS* ,// ,4X, 
-• I • , 7X, *Z(  I  )  ',//) 

UR  I TE ( 6  » 50  l  )  (  I, E(I), 1=1, NS2) 

501  FORMAT (  IX,  I4,E16.8) 

ZAP!  I)  ARE  THE  NON-DIMENSIONAL  STRENGTHS  OF  THE  SOURCE-SINK 
PAIRS  IN  BODY  A  THAT  ARE  SOLVED  FOR  IN  SUBROUTINE  SIMQ. 

ZBP(I)  ARE  THE  NON-DIMENSIONAL  STRENGTHS  OF  THE  SOURCE-SINK 
PAIRS  IN  BUOY  B  THAT  ARE  SOLVED  FOR  IN  SUBROUTINE  SIMQ. 

EVALUATE  CP  AND  BC  ALONG  BODY  A 

00  303  J=  1  , N S 
UX ! J ) =0 .0 
UR! J )=0.0 
BC< J)=0.0 
CP! J )=0.0 
X ! J ) =XA ( J ) 

R ( J ) =R A ( J  ) 

303  CONTINUE 
C 

DO  4400  J -  1 , NS 
DO  405  1  =  1, NS 

PA1=  SORT!  (  X ( J )-XA  (  I ) )**2+R!J)**2) 

PA  2=  SORT! !X!J)-XA!  I) )**2+!DELA( I  )-R!J) )**2 ) 

PB I =  SORT!  !  X ( J )-XB ! I ) ) **2  + I R { J ) -YO ) **2 ) 

PR2=  SORT!  ! X ! J ) -XB (  I  )  ) **2  +  l R ( J ) ♦DELB!  I  )-Y0)**2 ) 

C 

AX= ( X ( J )— XA ( I)  )/HAl**3 

BX  =  ( X( J ) -X A ( I ) )/PA2**3-!X( J)-XA ( I ) ) /PA  1**3 
C  X  = ( X<J  )  —  X  B ( I ) ) /PE1**3 

DX= ( X ( J ) -X B !  I)  )/PB2**3-!X! J)-XB< I  )  )/PBl**3 

C 

AR=R ! J ) /PA  1**3 

BR  =  (R(J)-DELA(I ) ) /P A2* * 3-R I J ) /P A 1 ** 3 
CR=(R(J  )-Y0)/PBl**3 
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ooo  rv>  o  O  o  o  o 


DR=(RtJ)+DELHCI ) -Y0 ) / PB2* *3- ( R ( J ) - VO > /PB 1**3 


UX(J)*UX(J)47A(I)*AX+ZAP(  I  )*BX+ZB(  I  )*CX+ZBP(  l)*DX 
UR<J)*UR( J  )  +  ZA<  I  )*AR  +  ZAPI I) *BR+ZB (II *CR+ZBP( I ) *0R 
405  CONTINUE 

BC(J)=UR(J)/(1.G+UX(J)  ) 

CP( J )=-2.0*UX< J)-UX< J)**2-UR(J)**2 

4400  CONTINUE 

WR  I  TE ( 6  » 60C )  XO»YO 

600  FORMAT (  1H1 ,  •  BODY  A  XO  *  '.F10.5,'  YO  =  SF10.5/ 

-•I  XA  RA  DRDXA  BC 

-UX  UR  CP  Z A  ZAP'//) 

WRITE (6, 700) ( I,XA( I ) , R A (I) , DRDXA ( I ) , BC ( I ) , UX I  I ) ♦ UR ( I ) 

-  »  CP(I),  Z  A  (  l  ) ,  ZAP! I ) ,  I  *  I *NS  ) 

700  FORMAT! IX, I2,9F13.5/) 

EVALUATE  CP  AND  BC  ALONG  BOOY  B 

DO  304  J= 1 ,NS 
UXI J)=0.0 
UR ( J ) -0 • 0 
BC ( J ) -0 • 0 
CP ( J ) =0 .0 
X! J)=XB(J) 

R ( J ) =YG-RB ( J ) 

304  CONTINUE 
C 

DO  5500  J=1,NS 
DO  505  1=1, NS 

P A 1=  SORT ( (X(J >-XA(I) )**2+R( J)**2> 

PA2=  SORT ( (X(J)-XA(I) )**2+(DELA(I)-R( J) )**2) 

PB  1=  SQRT  (  (X(J)-XB(I))**2MR(J>-Y0)**2) 

PB2=  SURT ( { X  C  J )-XB( I ) ) **2+( R ( J ) ♦DELB (  I  )-Y0)**2) 

C 

AX  =  (X(J )  —  X  A (  I  ) )/PAl**3 

BX  =  I X( J )-XA ( I ) )/PA2**3-(X( J)-XA( I )) /PA  1**3 
CX= I  X ( J ) -xe (  I)  )/PBI**3 

OX*(X(J)-XB(  I) )/PB2**3-(X( J)-XB( I ) )/PBl**3 
C 

AR=R ( J ) /PA  1**3 

BR  =  (R( J)-DELA( I ) ) /PA2** 3-R ( J ) /PA  I** 3 
CR=(K{ J )-YC)/PBl**3 

DR=( R ( J J+DELB! I )-Y0 ) /PB2**3- ( R { J )- YO ) /PB 1**3 
C 
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r>  o  o  o 


ux ( J )=UX( J )+ZA( I)*AX+ZAP( I i*BX+ZB(I )*CX+ZBP( I)*DX 
UR ( J )=UR( J  )  +  ZA  {  I )*AR*ZAP( I ) *BR+ZB ( I )*CR+ZBP( I  )*DR 

C 

50 5  CONTINUt 
C 

BC(J)=URI J  > / ( 1 . 0  +  UX  C  J  >  ) 

CP( J )=-2.C*UX( J)-UX( J)**2-UR( J)**2 

C 

*3500  CONTINUE 

WR I TE ( 6  *  8CC )  XO,YO 

RCO  EURMATUH1,'  BOOY  B  XO  =  *,F10.5,»  YO  *  *,F10.5/ 

-•  I  XB  RB  ORDXB  BC 

-UX  UR  CP  ZB  ZBP ' // ) 

WRITE  (6, 700)  U,XB(  I),RB(I  )  »  ORDXB  (  I  )  ,BC  ( I  )  ,UX  ( l  )  ,UR  (I) 

-  ,  CP(I  )  ,  ZB ( I ) »  ZBP  (  I  )  t  1=1, NS  ) 


EVALUATION  OF  UX,  UR,  UT ,  CP,  AND  BC  AROUND  BODY  A 
T  =  90  •  0 

24  00  20  J  =  1 »  NS 
UX(J)=0.00 
UR (J )=U.00 
U  T ( J ) =0 . 00 
BC I J ) *0.00 
CP(J)=0.00 
X ( J ) =XA ( J  ) 

R ( J ) =RA ( J ) 

20  CONTINUE 
C 

DO  21  J= l , NS 
DO  22  1=1, NS 
C 

PXA=  X( J)-XA ( I ) 

PXB=  X(J)-XBU) 

PZ=  COSO ( T  )  *R  <  J ) 

C 

P A 1  =  SORT (  PXA**2+PZ**2+(  SIND(T)*R(J)  )**2  ) 

PA2  =  SORT!  PXA**2'*'PZ*  +  2+  I  S  INDIT  )*R{  J  )-DELA(  I )  )**2  ) 

PR  1  =  SCRT (  PXB**2+PZ**2d-SIND(T)*R(  J  )  +Y0)**2  ) 

PB2  =  SCR  T (  PXB**2+PZ**2+ (-S I ND  < T )*R( J ) -DELB ( I )  +  Y0) **2  ) 

C 

AX  =  PVA/PA  1**3 
RX=PX4/PA2**3  -  AX 
CX=PXB/PBl**3 
0X=PXH/PB2**3  -  CX 
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AR  =  R(J)/PA1**3 

HR  =  <K( J)-SIND(T)*DELA( I ) )/PA2**3  -  AR 

CR=  - ( S IND ( T ) *YO-R (J) ) /PB1**3 

DR=-(SIND(T)*{Y0-DEL6(I))-R(J))/PB2**3  -  CR 

C 

AT  =  0.0 

BT  =  -CnSO(T)*  DELA! I )  /PA2**3  -  AT 

CT  =  -COSO(T)*  YO  /PB1**3 

DT  *  -CCSC ( T  )  * ( YO-DELB ( I))/PB2**3  -  CT 

C 
C 

UX(J)=UX(J)  ♦  Z  A  ( I )  *AX  ♦  ZAP  (  I  )  *BX  ♦  ZBin*CX  ♦  ZBP(I)*DX 

UR ( J ) =UR ( J )  +  ZA ( I ) *AR  +  Z AP { I ) *BR  ♦  ZB(I)*CR  ♦  ZBP(I)*DR 

22  UT  (  J  )  SUT  (  J  )  ♦  Z  A  (  I )  *AT  ♦  ZAPUJ^BT  ♦  ZB(I)*CT  ♦  ZBP(I)*DT 

C 

c 

BC ( J ) -UR ( J ) / ( I.O+UX(J) ) 

21  CP ( J ) =-2.0*UXl J ) -UX! J )**2-UR I J )**2-UT ( J )**2 

C 

c 

THE  T  A=T 


WRITE (6,25 ) 

XO,  YO, 

THETA 

FORMAT! 1H1, 

•  80DY  A 

XO  *  * ,F5. 3 , ' 

YO  = 

• ,F8.5,» 

THETA  *•, 

F6. 2/ 

•  I 

XA 

RA 

ORDXA 

BC 

UX 

UR 

CP 

ZA 

ZAP 

UT'//) 

WR I TE ( 6 , 333  ) (I  ,  XA (I > ,RA < I ) , DRDXAII) , BC (I)  ,UX < I ) , UR (  I ) 
-  ,  CPU),  ZA(I),  ZAP  (  1  )  ,  UT  (  I  )  ,  Is  1 , NS  ) 

333  FORMAT! IX, 12, 10F13.5/) 

IF(T .LE.-90.0)  GO  TO  23 
T  *  T  -  30.0 
GO  TO  24 
23  CONTINUE 
RETURN 
END 
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SGBROUT I Nfc  GEOM(RAD,X,OXfNS,SLOPE ) 

GEOMETRY  FOR  M117  BOMB 

REAL  LI ,L2,L3,L4,L5 
INTEGER  FINS 

DIMENSION  RADI  31 ) , X ( 3  l  ) , SLOPE  I  3 1 ) 

L  1  =  1 . 32288 
L2=2. 54375 
L  1=4 . 35G00 
L  4  =  4 . 9*t000 
L  5  =  6 . 30 
R  ADT  =  0. 262 
T AM  =0.  7559289 
TAN2=0. 1 316500 
NSPl=NS+l 
X ( 1 >  =DX 
F  I  NS  =2 

F I  NS  =  1  CORRESPONDS  TO  SIMULATED  FIN  BODY 
F 1 NS=2  CORRESPONDS  TO  STRAIGHT  TAPER  TAIL  BODY 

DO  9  I  =  2 » NSP  1 
X( I )=X{ I-i l+DX 

NOSE  SECTION 

IF(XII-l).GT.Ll)  GO  TO  l 

RADI I-l >=  SCRTC4.0-IXI I-l 1-1.32 )**2) -1.500 
SLOPE! I - 1 ) =- (XI I-l)- 1.32)/  SORT  1 4 1  X ( I- 1 )- 1 . 32 ) **2 ) 

go  rn  g 

MID-SECTION 

1  I  F  I  X  I  I-D.GT.L2)  GO  TO  2 
KADI  I-l  >=0.50 
SLOPE  I I-i )=0.0 
GO  TO  9 

STRAIGHT  REGION  OF  TAIL 

2  I F  I X  I  I-D.GT.L3)  GO  TO  3 
KADI  I-l  )  =0.5-1  XI  I-D-l  2)*TAN2 
SL0PE(I-l)=-0. 13165 
GO  TO  9 

3  GO  TO  (4,7) , FINS 
C 

C  SIMULATEC  FIN  BODY 
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4  IMXU-1J.GT.L4)  GO  TO  5 
RADI  1-1  )=RADT+(Xl  I - l)-L 3 ) * < C. 240/ ( L4-L 3 )  ) 

SLOPE ( I - L ) =C .240/ { L4-L  3 ) 

GO  TO  9 

5  IF  ( X  C  I-l»  .GT.L5)  GO  TO  6 
RAO< 1-1 )=0.502 
SLOP  £ ( 1-11=0.0 
GO  TO  9 

6  RAD(  1-11=0.0 
GO  TO  9 

STRAIGHT  TAPtR  TAIL 

7  I F ( X ( I-1J.GT.L5)  GO  TO  8 
RAO (  I-l  >=0.5-TAN2*(XU-l)-L2) 

SLOPE! l-l )=-0. 13165 

GO  TO  9  ** 

8  RAD(  1-1 )=0.0 
SLOPE! 1-1 ) =0.0 

C 

9  CONTINUt 
C 

GO  TO  ( 10,  12), FINS 

10  WKITE(6,ll) 

11  FORMAT ( 1H1 ,’ GEOMETRY  FOR  SIMULATED  FIN  BODY*//*  I 

-RAC(I)  SLOPE { I ) * / 1 

GO  TO  14 

12  WR I T  E ( 6 , 1 3  ) 

13  FORMAT ( 1H1 , ’GEOMETRY  FOR  STRAIGHT  TAPER  TAIL*//' 

-RAO(I)  SLOPE ( I  )  ’ / 1 

14  WRITEI6.15)  (I,X(I),RAO(n,SLOPEm,I  =  l,NS) 

15  FORMAT! 15, 3F10.5) 

C 

RETURN 

END 
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SUBROUTINE  S I  MO ( A,B,N,KS) 


DOUBLE  PRECISION  SIMQ 

DIMENSION  All) ,B(1) 

DIMENSION  A ( 1 ) »  B ( 1 ) 

FORWARD  SOLUTION 

TOL-O. 0 
KS  =  0 
JJ=-N 

DO  65  J  =  1 »  N 
J Y= J+ 1 
J J=JJ*N+1 
B I GA=0 
I T=J J- J 
DO  30  I  =  J  » N 

SEARCH  FOR  MAXIMUM  COEFFICIENT  IN  COLUMN 


I J- I T+ I 

I F (  ABSIBICA)-  ABS( A( I J ) ) )  20,30,30 
20  B  IGA=A ( IJ ) 

I  MAX  =  I 
30  CONTINUE 

TEST  FOR  PIVOT  LESS  THAN  TOLERANCE  (SINGULAR  MATRIX) 

I E (  AHS(BIGA)-TOL)  35,35,40 
35  KS  =  1 
RETURN 

INTERCHANGE  RCWS  IF  NECESSARY 

40  I  l  =J  +N*  (  J-  2  ) 

I T=I MAX- J 
DO  50  NK  =  J  ,N 
I  l  =  I 1  +  N 
12=1 l+IT 
SAVE=A(  1 1  ) 

A(  U)=A(  12) 

A ( 12 ) =S AVE 

DIVIDE  EOUATICN  BY  LEADING  COEFFICIENT 


50  All l )  =  A{ I  1 )/BIGA 
SAVE=B( IMAX) 
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HI  IMAX)=B(  J) 

B ( J ) =SA VE /B IGA 

ELIMINATE  NEXT  VARIABLE 

IF(J-N)55,70,55 
55  igS=N»(J-l) 

DO  65  I X* J Y « N 
IXJ=IQS+IX 
I T=J-I X 
00  60  JX=JY,N 
IXJX=N*( JX-I ) ♦ IX 
J JX= I X JX+ 1  T 

60  A ( I  X J  X ) =A { IXJX)-(A(  IXJ)*A( JJX) ) 
65  B( IX ) =B ( IX)-(B(J)*AI I X J ) ) 

BACK  SOLUTION 

70  NY =N- 1 
I TSN*N 

DO  80  J  *  1 »  N Y 
I A=I T- J 
I B=N- J 
I  C=N 

DO  80  NK  =  I »  J 
B(IB)=B(IB)-A(IA)*B(IC) 

I A=I A-N 
80  IC=IC-1 
RETURN 
END 
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19  ABSTRACT 

The  mutual  aerodynamic  interference  problem  for  two  axisymmetric  bodies  has  been 
analyzed  using  the  image  system  technique.  In  order  to  facilitate  this  analysis, 
it  has  been  assumed  that  small  perturbation  solutions  are  valid.  It  is  further 
assumed  that  the  external  stores  are  slender  bodies  and  that  the  mutual  inter¬ 
ference  can  be  analyzed  by  first  assuming  a  cross-flow  solution.  The  image 
system  in  the  cross-flow  plane  consists  of  source-sink  pairs  appropriately  located 
by  using  the  Milne-Thomson  circle  theorem.  The  actual  three-dimensional  source- 
sink  pairs  are  displaced  from  the  body  axis  according  to  the  cross-flow  image 
system.  Their  strengths  are  then  determined  by  the  Rankine  method.  Good  agree¬ 
ment  has  been  found  between  the  theoretical  and  experimental  results. 
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